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Abstract.  Our  interest  in  the  stability  analysis  of  the  high-order  continuum  (HOC)  equations  is  motivated  by  the 
relevance  to  the  development  of  a  hybrid  method  combining  such  equations  with  the  Direct  Simulation  Monte-Carlo 
(DSMC)  technique  for  the  computation  of  hypersonic  flows  in  all  regimes  -  continuum,  transition,  and  rarefied.  The 
hybrid  approach  allows  the  effects  of  thermophysics  (thermal  and  chemical  non-equilibrium)  and  turbulence  to  be 
included  much  more  easily  than  in  other  approaches,  and  can  easily  be  developed  into  a  robust  and  efficient 
engineering  tool  for  practical  3D  hypersonic  computations.  Stability  characteristics  of  model  HOC  equations  when 
subjected  to  small  disturbances  are  investigated.  We  explore  the  feasibility  of  simplified,  yet  accurate  and  numerically 
stable,  versions  of  the  HOC  equations  and  extend  our  previous  work2  to  include  multidimensional  Burnett  equations, 
with  the  specific  example  of  the  Augmented  Burnett  models.  The  latter  is  shown  to  have  a  much  wider  stability  regime 
than  Lumpkin’s  model. 


INTRODUCTION 

Hypersonic  flows  about  space  vehicles  produce  flow  fields  with  local  Knudsen  numbers,  Kn ,  which  may 
lie  in  all  the  three  regimes  -  continuum,  transition,  and  rarefied.  The  Navier-Stokes  (NS)  equations  and  the 
direct  simulation  Monte-Carlo  (DSMC)  methods  can  accurately  and  efficiently  model  flows  in  the 
continuum  and  rarefied  regimes,  respectively.  The  kinetic  approach  considers  an  ensemble  of  small 
particles  or  molecules  whose  distribution  function  can  be  determined  as  a  solution  of  the  Boltzmann 
equation,  while  the  continuum  approach  is  based  on  the  representation  of  the  gas  as  a  fluid  continuum 
governed  by  the  mass,  momentum,  and  energy  conservation  laws.  Though,  theoretically,  the  kinetic 
approach  is  appropriate  for  simulating  gas  flows  in  any  regime,  in  practice,  it  can  require  large  computer 
resources  if  the  gas  flow  is  dense.  DSMC  remains  the  most  efficient  numerical  technique  for  solving  the 
Boltzmann  equation  [1].  It  enables  the  computation  of  flows  with  high  Knudsen  number.  Nevertheless, 
DSMC  computations  are  still  too  expensive  in  many  cases,  especially  for  3D  engineering  applications. 
Although  a  rather  efficient  tool  for  supersonic  and  particularly  hypersonic  flows,  the  DSMC  procedure 
becomes  more  resource-consuming  for  low  Mach  number  subsonic  flows,  due  to  difficulties  with  boundary 
condition  implementation  on  subsonic  inflow/outflow  boundaries.  Furthermore,  obtaining  gas  interactions 
with  DSMC  is  a  difficult  task.  The  continuum  approach  is  much  cheaper  and  more  versatile  in  these 
regards.  There  is,  therefore,  a  strong  motivation  for  its  utilization  at  the  low  Kn  values.  The  traditional 
continuum  model  is  based  on  the  Navier-Stokes  equations,  which  are  the  first  order  approximations  to  the 
Boltzmann  equation  with  respect  to  the  (small)  parameter  Kn.  Coupled  with  no  velocity  slip/no  temperature 
jump  solid  wall  boundary  conditions,  they  are  valid  if  the  Knudsen  number  is  small,  say,  less  than  0.001. 
More  rarefied  flows  should  be  described  using  the  Navier-Stokes  equations  with  velocity  slip/temperature 
jump  boundary  conditions. 

The  flows  in  the  transitional  regime  require  higher-order  continuum  (HOC)  models;  the  most  well-known 
being  the  Burnett  equations,  obtained  as  second  order  approximations.  Though  there  are  some  difficulties 
with  the  stability  of  their  solutions  and  the  development  of  relevant  solid  wall  boundary  conditions,  recent 
advancements  [2]  allow  the  consideration  of  the  (properly  modified)  Burnett  equations  as  a  potential 
continuum  model  for  transitional  flows.  In  recent  years,  Burnett  equations  have  been  successfully 
employed  to  compute  3D  hypersonic  flows  in  continuum- transition  regime  [3],  although  it  has  been 
difficult  to  compute  flows  for  Kn>  1. 

The  other  high-order  continuum  (HOC)  equations,  such  as  Eu’s  [4]  and  Grad’s  13-moment  equations  [5], 
are  significantly  more  expensive  to  compute  than  the  Burnett  equations,  and  have  been  tested  only  for  ID 
and  for  2D  geometrically  simple  problems.  Another  approach  is  due  to  Aristov  and  Tcheremissin  [6], 
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wherein  the  left-hand  side  of  the  Boltzmann  equation  is  solved  with  a  finite-difference  scheme  and  a 
special  quadrature  formula  is  employed  for  the  collision  integral  on  the  right-hand  side.  This  method  has 
recently  been  applied  to  solve  2D  problems  involving  a  mono-atomic  gas.  Application  of  the  approach  to 
gases  with  internal  degrees  of  freedom  is  problematic  at  the  moment  because  of  the  difficulty  with  the 
inclusion  of  chemical  reactions.  We  therefore  investigate  the  Burnett  equations  for  use  as  the  HOC 
component  of  our  hybrid  procedure. 

In  this  paper,  we  examine  the  stability  of  a  few  versions  of  the  Burnett  equations,  to  guide  the  selection  of 
the  model  for  a  robust  hybrid  procedure  for  hypersonic  flows.  In  specifics,  we  will  examine  Lumpkin’s 
simplified  model  [7],  the  generalized  Burnett  equations,  and  the  Augmented  Burnett  equations  for  their 
relative  stability  characteristics.  Although  the  stability  analysis  of  the  Burnett  equations  have  been  reported 
in  the  literature,8  the  studies  did  not  include  Lumpkin’s  model  or  considered  the  presence  of  rotational 
temperature  in  the  analysis. 


ANALYSIS  OF  LUMPKIN’S  SIMPLIFIED  BURNETT  EQUATIONS 


The  one-dimensional  equations  are  considered  for  the  analysis  of  Lumpkin’s  simplified  model,  with 
TU  =  8  and  q_(q  +  2  ^  +  2  ^  3  in  the  standard  Burnett  equations.2,3,8  The  resulting  equations  can  be 
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written  as  follows: 
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and  E  =  I +  2rt  +  u2\  P  =  PRT  Tt ,  jU  ,  and^  are  the  total  energy/unit  mass,  hydrodynamic 

pressure,  rotational  temperature,  translational  temperature,  molecular  viscosity,  and  the  rotational  collision 
number.  Note  that  T  R  is  the  relaxation  time  for  rotational  energy,  so  that^  T  =  z  \  where  T  c  is  the 

* c  RW) 

mean  collision  time. 


Linearization  of  the  Equations 

Consider  a  diatomic  gas  in  equilibrium  with  density  p0 ,  pressure  p0 ,  translational  temperature  T;()  and 
rotational  temperature  TR 0 .  The  gas  is  subjected  to  small  perturbations  defined  as  the  non-dimensional 
variables: 

<=P-Po  ,  r,=  Tt-Tt0  ,r  ,=  Tr  ~Trq  ,  M,=  u  ,t'= L - 

Po  (  To  Tr  0  M0/p0  L0 

L  _  Po  .  Note  that  TR0  =  Tt0  is  assumed. 

Po  y/RTto 

The  linearized  equations  can  be  written  as: 
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and  j[  =  [0  0 --C 5Cf  >  where  “T”  denotes  the  transpose  operator. 

Solution  of  the  Linearized  Equations 

Consider  the  homogeneous  portion  of  the  linearized  equations  and  assume  V'=  Vel,x' e*' ,  where 
<p  =  a  +  ip  and  co  _  2;r  .  Since L  =  fi0  =  Q  ?g3;l ,  we  have  a  =  4  92£ = 4  92K  ,  where  Kn  is  the 
LIL°  °  P0W0  '  '  L  ~  ' 

Knudsen  number.  The  dispersion  relation  is 
det|^7  +  icoLl  -  cd2L2  +  L0 1  =  0 

which  can  be  simplified  to 
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This  equation  was  solved  using  MATLAB  to  determine  the  stability  boundaries  for  z^=  4,  10,  18,  and  23. 
Note  that  Lumpkin  recommended  \%<zR  <23  for  his  model  and  that  Jean’s  equation  has  been  used  in  the 
source  term. 

The  stability  boundaries  are  shown  in  figures  1  and  2  for  various  values  of  zR  •  Regions  with  a  <  0  (on 
the  x-axis)  are  stable,  whereas  regions  with  a  >  0  are  unstable,  y  is  taken  as  1 .4  for  a  diatomic  gas. 


FIGURE  1.  Stability  boundaries  of  Lumpkin’s  simplified  Burnett  model  for  Z R  =4  and  10 
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FIGURE  2.  Stability  boundaries  of  Lumpkin’s  simplified  Burnett  model  for  zR=  18  and  23 


ANALYSIS  OF  THE  GENERALIZED  3D  BURNETT  EQUATIONS 


The  governing  equations  considered  for  the  stability  analysis  of  the  3D  Burnett  equations  with 
translational  and  rotational  thermal  non-equilibrium  are  as  follows: 


,  djpu)  t  d(pv)  ,  djpw)  Q 
dt  dx  dy  dz 


!  d  (j 
dt  dx 

d{pv)  +  Z( 

dt  etc 
o(pw)_+d_l 
dt  dx 


dy 


dr  drn,  dr 


■f( 

dy 

f( 

dy 


i1 


-( 

dz 


dx 

dy  dz 

,)  =  _£5s 

,  8tw  dryz 

dx 

dy  dz 

dx 

dy  dz 

^^■  +  -^-{pu  +  pEu)  +  -^-{pv  +  pEv)  +  -^-{pw  +  pEw)  = 
dt  dx  dy  dz 


iptu  i-  (  l'|+  '  '  {pe,  u  |=  - 

dt  dx  dy  dz 


ut  +  vt  +wt  )h - fwr  +  vr  +wr 

‘  xy  yy  yz)  Qz  V  xz  yz  rv  y  2 


f 


(12) 

(13) 

(14) 

(15) 


t  dZl  !  dg; 
dx  dy  dz 


dqJL+dqJL+dqL 

dx  dy  dz 


P^r 

Zr 


(16) 


Note  that  the  last  equation  is  for  the  rotational  energy  and  that 
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=  —(6T  -T  ),  _^/is  mean  collision  time,  as  before,  so  that 
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The  Linearized  Equations 

Using  a  similar  non-dimensionalization  scheme  as  above,  the  equations  can  be  written  as 
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A. 

B. 

C. 


L3,  M3,  N3j  M4  =  0  ->  Navier-Stokes  equations  with  translational  and  rotational  non-equilibrium. 
This  limit  has  been  analyzed  for  various  values  of  ZR,  and  they  are  known  to  be  stable  for  ZR  =  0. 
M4  =  0  ->  conventional  Burnett  equations  with  translational  and  rotational  non-equilibrium.  The 
equations  are  known  to  be  unstable  for  all  ZR  including  ZR  =  0. 

M4  ^  0  ->  Augmented  Burnett  equations.  We  studied  their  stability  for  various  values  of  ZR.  For 
ZR  =  0,  they  are  known  to  be  stable  if 


(23) 


We  will  consider  0  <  ZR<  23  in  Case  C.  By  changing  the  values  of  oo7,  06,  and  07,  we  will  try  to  extend 
the  stability  for  the  largest  ZR  value. 
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FIGURE  3.  Stability  boundaries  of  the  Augmented  Burnett  equations  for  ZR  =4  and  23 


RESULTS  SUMMARY 

It  is  apparent  that  the  Lumpkin’s  equations  are  unstable  to  small  perturbations  in  a  quiescent  fluid  when 
ZR  >0  and  stable  otherwise.  Although  the  simplified  model  seems  to  work  well  in  some  cases,2  it  will  be 
necessary  to  use  either  the  Augmented  or  BGK-Burnett  model  to  include  the  rotational  non-equilibrium. 
However,  detailed  stability  characteristics  of  these  equations  are  of  interest  and  will  be  carried  out  in 
further  studies.  Note  that  ^  =  4.92^ ,  and  that  the  stability  boundaries  have  been  determined  by  varying  Kn 

from  0  to  1  using  100  points.  We  have  also  examined  the  stability  of  the  ID  Augmented  Burnett  equations 
for  ZR  =4,  18,  and  23.  The  results  for  ZR  =4  and  3  are  shown  in  figure  3.  The  results  for  the  Case  ZR=  18  are 
similar  to  those  for  ZR  =  23  and  are  therefore  not  shown  in  this  paper.  As  the  figures  show,  we  have  found 
that,  with  the  appropriate  coefficients,  the  Augmented  equations  are  stable  for  ZR  values  up  to  23.  This 
suggests  the  superiority  of  the  Augmented  Burnett  equations  over  Lumpkin’s  simplified  model.  Hence,  the 
former  might  be  more  appropriate  for  subsequent  work  on  hybrid  HOC/DMSC  procedures.  We  have 
incorporated  a  few  physically  realizable  and  computationally  stable  versions  of  the  Burnett  equations  into 
LAURA  [9]  and  combined  this  with  a  DSMC  procedure.  Details  are  available  in  [10]. 
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